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We have analysed a number of experimental geometries (two and three-beam geometry, general 
systematic row setup, zone axis orientation) for the electron magnetic circular dichroism measure- 
ments (EMCD). The effect of asymmetries in those setups is described in detail, explaining the 
advantages of the three-beam geometry. We observe a complicated thickness and orientation depen- 
dence of the dichroic signal and asymmetries. Therefore we stress the importance of the simulations, 
which should accompany the EMCD experiments. 
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I. INTRODUCTION 

The recent discovery of the possibility to observe the 
magnetic circular dichroism in the transmission electron 
microscop^ opened a new route to atom-specific mag- 
netic studies with an unprecedented spatial resolution. 
This technique was named electron magnetic circular 
dichroism (EMCD). Although principle of this technique 
shares many common points with x-ray magnetic circular 
dichroism (XMCD), dynamical electron diffraction has a 
pronounced effect on the observed spectra. This brings 
many subtleties for the interpretation of experiments and 
the extraction of physical quantities. In order to turn the 
technique into a quantitative probe of magnetic proper- 
ties, it is thus necessary to analyse possible sources of in- 
accuracies connected with dynamical diffraction. We will 
deal with so far the only viable EMCD approach, in which 
the crystal itself is used as beamsplitter^ and the EMCD 
signal is obtained as a difference of spectra measured at 
two detector positions. We concentrate here on the role 
of symmetry of the experimental setup geometry, i.e., rel- 
ative orientation of incoming beam vs sample orientation 
vs detector positions. Note that yet before the EMCD 
was discovered, signals of odd and even symmetries were 
separated in an electron energy loss experimeniP, though 
relation between the signal of odd symmetry and mag- 
netism were not yet noticed. 

In a recent publication'^ we have described simulations 
of electron energy-loss near edge structures ( ELNES ). 
This theory is based on the Bloch wave approacfP^I^ to 
the dynamical diffraction effects, including higher-order 
Laue zones. Matrix elements of inelastic transitions are 
evaluated from first- principles electronic structure. This 
method was applied to study EMCD in the experimental 
geometry proposed in the original papei^. Those simu- 
lations were done within the systematic row approxima- 
tion. Inclusion of high order Laue zones provides more 
complete, richer description leading to qualitatively new 
observations presented in this work. 

An important theoretical development for quantitative 



use of th e EMCD technique was the derivation of sum 
rule^^Hni which relate the measured spectra to spin (ms) 
and orbital (m/,) magnetic moments of a studied atom in 
the sample. Among the first quantitative measurements 
of the mi/mg ratio are Refs. lO'll' and, particularly, 
Ref . 1121 where a method of data treatment and statistical 
analysis was described. 

The sum rules require two measurements in symmet- 
rical conditions for determination of the orbital to spin 
moments ratio. When the symmetry condition is ful- 
filled, sum rules revealed a flexibility of the EMCD tech- 
nique allowing to measure the dichroic signal in the whole 
diffraction plane. Consequences of these findings are the 
main points of this paper. 

After describing computational methods (Sec. [ll]) we 
study in detail the originally proposed experimental ge- 
ometry (Sec. fill - the two-beam case (2BC). Maps of the 
dichroic signal are presented. We demonstrate that 2BC 
actually does not provide two strictly symmetric mea- 
surements and analyse the error introduced to the ex- 
tracted TOi/ms ratios. An improved experimental setup 
based on three-beam case (SBC) geometry is discussed in 
Sec. |fV[ In Sec. |V] we study the general systematic row 
geometry. Finally, in Sec. |Vf| we analyze the zone axis 
(ZA) orientation and its suitability for EMCD measure- 
ments. 



II. METHOD OF CALCULATIONS 

In Ref. |4]we have described our theoretical approach to 
simulate general orientation-sensitive ELNES experiment 
from first principles. These methods will be fully utilized 
in this work. Here we present calculations for bulk bcc 
Fe including higher order Laue zones, Weickenmeier-Kohl 
scattering potentiald^ and Debye-Waller factor O.SSA^ 
for Fe at room temperatur^^^. The energy of the elec- 
tron beam was set to SOOkeV and the sample thick- 
ness varied between 5 nm and 50 nm. For solution of 
the eigenproblem for both incoming and outgoing beam 
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FIG. 1: (color online) The Thales circle circumscribed to ex- 
citations and G is a set of detector positions, for which 
the two momentum transfer vectors q and q' are perpendicu- 
lar. The detector positions denoted as (A) and (B) fulfill also 
I q| = I q' I . These two detector orientations are usually quoted 
as the Thales circle positions. 



we use basis consisting of about 700 beams (maximum 
Wg = ^gSg = 10^, where and Sg are the extinction 
distance and excitation error, respectively). 

Compared to our previous work, Ref. El we have sig- 
nificantly improved the summation over the Bloch co- 
efficients. Details of implementation and study of the 
convergence are given in Appendix [X] here we list only 
the convergence criteria. A dipole-type MDFF asymp- 
totics is assumed in estimating the term magnitudes and 
all terms exceeding 10~^ are included in the summation. 
For zone axis simulations an increased limit 3 x lO^"' was 
used to reduce computational costs. 

For calculation of matrix elements S{q,q',E) - the 
mixed dynamic form factors (MDFF)- - we used density 
functional theory. Spin-orbit interaction was included 
and magnetic moments were aligned with incoming beam 
direction. Particularly, M||(016) to simulate ^10° tilt of 
the sample within systematic row setups with G = (200). 
In zone axis simulations, M||(001). We use Rayleigh 
expansion of the Coulomb operator to the first order 
(A = A' = 1 approximation, where A is the order of the 
terms included in expansion), which keeps the computa- 
tional efficiency of dipole approximation while producing 
much more accurate results for larger momentum trans- 
fer vectors. We refer the reader to the above-mentioned 
paper- for more details on the simulations. 



III. TWO-BEAM CASE 

We will describe in more detail the originally pro- 
posed geometry of the EMCD experiment, as applied in 
Refs. lll4l This geometry is based on exciting a systematic 
row of reflections (determined by the reciprocal lattice 
vector G) by tilting the incoming beam ki„ a few degrees 
from some highly symmetric direction, such as (001), 
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FIG. 2: (color online) Asymetry of the two-beam case. The 
incoming beam with wavevector kin is tilted with respect to 
the zone axis, which is perpendicular to the sample surfaces. 
The square grid is a sketch of the crystal lattice. Detected 
outgoing beams for two detector positions (A) and (B) are 
described by their respective wavevectors kout. The dotted 
light arrows are auxiliary vectors parallel with incoming and 
outgoing beam wavevectors used for construction of momen- 
tum transfer vectors q(A) and q(B) (thick arrows inside the 
crystalline lattice). 



within a plane perpendicular to the selected G . Then the 
incoming beam is tilted a few miliradians further along 
the systematic row in order to produce a two beam case 
(see Fig. 3 of Ref. H}. This is achieved by setting the 
Laue circle center of the incoming beam to Icc = G/2. 
We define momentum transfer vectors q = ko^t — ki„ and 
q' = q -I- G, with respect to the transmitted and Bragg 
scattered beams, respectively. By taking the excitations 
and G as the diameter of a Thales circle, then for all 
detector positions hout lying on this circle, the projec- 
tions of q and q' on the diffraction plane are perpen- 
dicular. A particularly symmetric situation is achieved, 
when |q| = |q'|, i.e., the detector placement in a direc- 
tion corresponding to the top (A) or bottom (B) point 
of the Thales circle, see Fig. [T| In a simple two-beam 
approximation, measuring at these detector positions is 
formally in the closest analogy to XMCD experiments.!^ 
The archetype EMCD experiment thus amounts to ac- 
quiring two spectra at these two detector positions. The 
difference of these spectra is attributed to the magnetic 
circular dichroism, because in the presence of inversion 
symmetry the A^IDFF changes sign of its imaginary part, 
S{q,q',E) = S{~q',-q,E) = S*{q,q',E),^ when mov- 
ing from detector position (A) to (B), see Fig. [l] The 
real part is assumed to stay unchanged and cancel out in 
the difference. Then the difference spectra would depend 
only on imaginary parts of contributing MDFFs. Energy 
integral of the imaginary part of MDFF was shown to be 
proportional to magnetic moments of the studied atoirP. 
Now we will demonstrate, that actually the (A) and 
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FIG. 3: (color online) Up-down relative diflerence maps of Fe- 
L3 edge calculated in 2BC as a function of thickness within 
(top row) systematic row approximation, and (bottom row) 
including HOLZ. The s-axis coincides with the systematic 
row and the origin of the coordinate system is at the position 
of the transmitted beam (in the middle). Tick labels of axes 
are multiples of the Bragg-reflection G'{2oo)- The same area 
of the diffraction plane is used in all subsequent figures. 

(B) detector positions are, strictly speaking, not sym- 
metric and therefore the signal measured at these two 
positions differs more than solely by reversing the sign of 
imaginary parts of contributing MDFFs. 

The assymetry itself appears due to the tilt of the crys- 
tal (usually approx. 5 to 10 degrees). The situation is 
schematically drawn in Fig. [2] Note that 1) the outgo- 
ing beam wave vectors kout for the two detector positions 
have different relative orienation with respect to the crys- 
tal lattice; 2) the same is valid for the momentum trans- 
fer vectors q(A)iQ(B) drawn inside the sample. These 
two facts are reflected in two sources of asymmetries: 1) 
different reciprocal wave Bloch fields at detector position 
(A) vs detector position (B), and, 2) difference in calcu- 
lated MDFFs, which goes beyond the mere change of the 
sign of the imaginary part. 

The difference in outgoing Bloch wave fields appears 
because of the different relative orientations of the Ewald 
spheres for observed outgoing beams with respect to the 
reciprocal lattice. Note though, that in the systematic 
row approximation (SRA), i.e., when we take into ac- 
count only Bragg reflections along the x-axis (perpen- 
dicular to the Fig. |2] plane), then the two eigensystems 
become equal. In such case, from the whole reciprocal 
lattice we consider only the systematic row and such 
simplifled geometry acquires a mirror symmetry plane 
defined by x and ki„. Nonetheless, this is an approx- 
imation and it is necessary to analyse its validity and 
appropriateness 

Strength of the latter effect depends on the anisotropy 
of the electron distribution surrounding the excited 
atoms. Note that this error also increases with larger 
energy losses, since larger energy loss shortens the kout 
more and the angle between related momentum transfer 
vectors increases. The effect itself introduces differences 
between the real parts of MDFFs, which were expected 
to be symmetry-related. As a consequence, the difference 
signal from positions (A) and (B) will contain not only 
the imaginary parts (the magnetic signal) but also the 
real parts of MDFFs. 



3 




-2-10 1 2-2-10 1 2-2-10 1 2-2-10 1 2-2-10 1 2 



FIG. 4: Apparent niL/ras ratio maps in the 2BC geometry as 
a function of thickness calculated within two models: system- 
atic row approximation (SRA, top row) and full calculation 
(bottom row). 



Figure [3] compares "up-down" difference maps calcu- 
lated within the SRA with calculation including higher 
order Laue zones (HOLZ calculation). Sec. |lT] They dis- 
play the antisymmetric part of the signal with respect 
to the systematic row axis - i.e., the difference of the 
diffraction pattern and its mirror image with respect to 
the systematic row. We will call them simply difference 
maps, contrary to names 'dichroic maps' or 'dichroic sig- 
nal' used in recent literature. The reason is that in the 
majority of the experimental geometries this difference 
signal is only an approximation to the magnetic com- 
ponent. Due to asymmetries these maps contain also a 
fraction of non-magnetic signal, which can introduce a 
systematic error into the mj^jms evaluation. Relative 
difference maps are difference maps nomalized by the 
symmetric part of the signal, i.e., the difference divided 
by the sum. 

There are several local maxima and minima in the 
diffraction plane. Fig. |3] Their positions and strengths 
depend on thickness and, in general, none of them is on 
the Thales circle. Remembering that the sum rules allow 
mj^/ms extraction at any place in diffraction plane, it is 
obvious that Thales circle positions are not the optimal 
detector position for measurement, when it comes to op- 
timizing signal-to-noise ratio (SNR). Actually, there is no 
general 'optimal detector position'. The optimal detector 
position is material, thickness and orientation dependent 
(see also Appendix [b]). This problem was studied in detail 
in Ref. TSifrom the SNR point of view. 

The main features of HOLZ calculation are reproduced 
by the SRA, see Fig. [3j However, there are some quali- 
tative changes in the HOLZ simulation, which are not 
present within SRA and the deviations increase with 
thickness. Two pronounced features appear in the HOLZ 
maps, which are not present in the SRA maps: Kikuchi 
lines and weakly excited Bragg spots away from the sys- 
tematic row. A Kikuchi line and its mirror image can be 
best seen at thicknesses 40 nm and 50 nm - the skew lines 
in the upper left and lower left corners of HOLZ maps. 
The weakly excited Bragg spots are almost invisible in 
diffraction patterns, but the relative difference maps en- 
hance their visibility. Presence of both these effects in 
relative maps is a sign of the asymmetry. 
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If these maps would originate solely from imaginary 
parts of MDFFs (a pure magnetic signal), it should be 
possible to obtain mi^/mg ratio at any position in the 
diffraction plane. This follows from the form of spin and 
orbital sum rules expressions, according to which the dy- 
namical diffraction effects cancel out in the m^/ms r atio. 
We tested this by evaluating the sum-rules expressiorP^^ 



ms 



2 A, 



(1) 



throughout the diffraction plane (here A^^ 3 is the differ- 
ence signal integrated over L2 or L3 edge, respectively) . 
On top of the usual sum rules approximations ,1^ by writ- 
ing this expression we assume that i) both spin and or- 
bital momentum vectors are parallel and ii) the magnetic 
dipole term is negligible compared to the spin moment. 

The results of these simulations are in the Fig.|4] where 
we compare such evaluation within SRA and HOLZ cal- 
culation. The "mi/ms ratio" is by no means constant 
in the diffraction plane, not even in SRA. Therefore, in 
the following text, we will refer to values in these maps 
as apparent rriL/ms ratio to distinguish them from the 
fnL/ms ratio obtained by electronic structure calcula- 
tions, which is 0.045 for the energy integration range used 
in our simulations (1.5 — 2.0eV above the Fermi level, see 
Appendix [C]) . 

There are several factors that affect the rriL/ms ratio: 
1) the assumption of equal radial wavefunctions for both 
2^1/2 E^iid 2p3/2 initial states made in derivation of sum 
rules, 2) slight difference in the direction of the spin and 
orbital moments due to magneto-crystalline anisotropy, 
3) magnetism of s-states and related p — > s transitions 
and p s/p ^ d cross-terms in inelastic transition ma- 
trix elements, and, 4) the asymmetry discussed above. 

The first effect can be easily quantified by performing a 
test calculation with equal initial-state radial wavefunc- 
tions for both edges. The fact that pi^2 and P3/2 radial 
wavefunctions are slightly different was observed to de- 
crease the apparent rriL/ms ratio by ~ 0.01 in the areas 
weakly affected by other sources of deviations. This ex- 
plains the mean value of 0.035 within regions close to the 
systematic row. 

The second effect depends on the strength of magneto- 
crystalline anisotropy, i.e., to what extent the (L) and (S) 
follow the external magnetic field and, particularly, what 
is the angle between (L) and (S). Considering the vector 
nature of sum rulea^, prefactor of the spin and orbital 
magnetic moments is actually a vector (p) and thus the 
ratio of type p • (L) /p • (S) may vary with varying p if 
(L) H (S). This is a tiny effect in bcc-Fe, which has ex- 
tremely low magnetocrystalline anisotropy. As a result, 
at M II (016) the angle between spin and orbital moments 
is below 0.1° and therefore the error thus introduced is of 
no practical importance. For systems with strong magne- 
tocrystalline anisotropy this effect might become visible. 
That might, in principle, allow to determine the angle 
between spin and orbital magnetic moment by following 



the changes of the m l /ms ratio throughout the diffrac- 
tion plane. However, detailed study of this effect goes 
beyond the scope of the present paper. 

The asymmetry effect is the most crucial one. It can 
lead to a quite large systematic error, when the mi,/ms 
value is extracted directly from 2BC maps further from 
the systematic row. 

The asymmetry within the SRA has a form of quasi- 
periodic pattern of vertical stripes. These modulations 
have a long period and low amplitude at low sample 
thicknesses, but as the thickness increases, the period 
shortens and amplitude increases. This effect can be 
traced doen to differences of elongations 7*^^^ of outgo- 
ing beam Bloch waves in upper and lower half-planes, 
which is caused by difference of projections of the outgo- 
ing wave vectors k(A) E^^id k^^) to the sample surface nor- 
mal (see also footnote in discussion of asymmetry under 
Fig. [2]). Thickness enters the differential scattering cross- 
section through factord^l e'*- which suggests that 
at larger thickness the scattering cross-section will be 
more sensitive to differences in elongations 7'-'-'. 

The asymmetry within the full HOLZ calculation is 
stronger and more complex. Notice, for example, that 
at the position of Kikuchi lines the apparent mj^jms 
ratio is estimated to be off by an order of magnitude. 
Deviations are stronger at detector positions, where the 
magnetic signal becomes weaker. Also, the deviations 
increase with increasing thickness, which puts an upper 
limit on sample thicknesses for EMCD measurements. To 
summarize, it is the asymmetry originating in the com- 
plete treatment of the dynamical diffraction, which has 
the strongest influence on observed 771^/7715 ratios in the 
2BC. A more detailed computational study of the asym- 
metry in the 2BC can be found in Ref. [T?l 

Still, there is a quite large area in the diffraction plane, 
where the 2BC geometry allows accurate 777^/7775 ratio 
measurement, Fig.|4] Though, due to dynamical diffrac- 
tion effects, the shape and position of this area is thick- 
ness dependent. Therefore it would be of advantage to 
accompany experiments in 2BC setup with simulations. 
Eventually, one can perform experiments at varying sam- 
ple thicknesses to identify regions with biased mi^jmg 
values. If a suitable thickness and detector position in 
diffraction plane is identified, this geometry might be pre- 
ferred for spectroscopic measurements of EMCD, because 
the two-beam orientation concentrates the magnetic sig- 
nal to the area inside and close to the Thales circle (com- 
pare with three-beam orientation. Sec. IV I. 



IV. THREE-BEAM CASE 

As we have shown, the two-beam case does not allow 
two independent strictly symmetrical measurements as 
required by the sum rules. An alternative experimental 
geometry, the three-beam case (3BC), is in principle free 
of these inaccuracies. Instead of exciting a 2BC using 
the incoming beam we set the Laue circle center of the 
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FIG. 5: (color online) Default 2BC (left) and proposed SBC 
(right) experimental geometries showing a tilted sample and 
detector positions (A) and (B). 

default alt. 1 alt. 2 
(1 00) (000) (000) 

deti^] (f fO) (f fO) (f fO) 

detif> (f fO) (f fO) (f fO) 

TABLE I: Parameters for 2BC and SBC geometries and 
Thales circle detector positions. All are given in terms of sys- 
tematic row vector G and Laue circle center (Ice) for incoming 
beam and detector position assuming the x axis parallel to G 
and z axis parallel to the zone axis. In SBC there are two 
possibilities for outgoing beams - 'above' (alt. 1) or 'below' 
(alt. 2) the systematic row, respectively. The right part of 
Fig. [5] shows alt. 1. 
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FIG. 6: (color online) Top row: maps of the L3 relative dif- 
ference signal of bcc-Fe in the SBC orientation constructed 
using the vertical mirror axis (solid line). Systematic row is 
shown as a dashed line and Bragg spots as small (blue) cir- 
cles. Bottom row: the rriL/ms ratio map constructed using 
the vertical mirror axis. 




FIG. 7: Map of the mL/ms ratio in the slightly misoriented 
SBC, Icc = (0.05,0,0) = 0.025G(2oo)- Top row: up-down 
difference (horizontal mirror axis), middle row: left-right dif- 
ference (vertical mirror axis), bottom row: double difference 
(both mirror axes). Values outside the range (—0.5,0.5) are 
white (if above 0.5) or black (if below -0.5), respectively. 



incoming beam to (000) . This geometry is conventionally 
called a three-beam case, because in the diffraction plane 
one sees a strong spot due to the transmitted beam and 
two G and — G Bragg spots of equal intensity. Other 
spots are typically much weak er. The SBC geometry was 
applied in recent experiments) "'^^ ! ^^ ^ 

If there exists a mirror symmetry plane normal to 
the systematic row (cf. mirror plane containing sys- 
tematic row in 2BC), this setup allows for two com- 
pletely symmetric measurements. Because the mirror 
symmetry operation M changes sign of pseudovectors, 
Mq X Mq' = — M(q x q'), the imaginary parts of MDFFs 
will change sign when going from (A) detector position to 
(B), see Fig. [s] The eigenvalue problems for the two de- 
tector positions are equivalent, because now the incoming 
beam remains within the mirror plane and the outgoing 
beams together with their associated Ewald spheres are 
mutually symmetric with respect to this mirror plane. 

Although the two-beam approximation and Thales cir- 
cle construction for the 2BC is an oversimplified model, it 
was very useful for getting an insight into the EMCD phe- 
nomenon and to develop formal analogies with XMCD 



techniqu^. It is thus instructive to show that a similar 
model can be constructed also for the SBC. In analogy 
with two-beam approximation for 2BC geometry, here 
we approximate the incoming beam by single plane wave 
(one-beam approximation) and set detector positions so 
that the outgoing Bloch field can be described by the two- 
beam approximation. It is possible to introduce similar 
Thales circle construction like in 2BC geometry, see fig.|5] 
and follow the analogy with XMCD in the same way as in 
the 2BC. The only difference is, that the 'mirror axis' is 
vertical and there are two Thales circles circumscribed to 
and G or and — G, respectively. Laue circle centers 
and Thales circle detector positions for both geometries 
are summarized in Table HI 

As mentioned above, the mirror symmetry plane cuts 
the diffraction plane in a direction perpendicular to the 
systematic row direction. Therefore to utilize the symme- 
try we evaluate difference maps of left and right diffrac- 
tion half-planes (cf. upper minus lower half-plane in 
2BC). If we would have perfect systematic SBC ori- 
entation, under dipole (or A = 1) approximation this 
difference map would contain solely a magnetic signal. 



This is apparent from Fig. 6] The rriL/ms rnaps are 
very smooth with only small variations of the apparent 
"^l/"^s ratio throughout the plane. These variations 
are caused pred ominantly by p — > d/p — > s transition 
cross-termg^nUl] Larger variations may appear at some 
particular thicknesses and positions - these are caused 
by reduced dichroic signal due to dynamical effects (typ- 
ically these are borderlines between areas of positive and 
negative dichroic signal). At these positions the m^/ms 
ratio is evaluated as a ratio of two very small numbers, 
i.e., becomes sensitive even to tiny deviations. 

In practical measurements, though, one never reaches 
a perfect SBC orientation. Either G or — G is slightly 
more excited than the other. Therefore it is instructive 
to test, how sensitive is this setup to small deviations. 
A combined theory-experiment study of this effect can 
be found in Ref. |22| Here we summarize our findings 
and present additional simulation results (see also Ap- 
pendix |d]) . 

By setting IcCm to (0.05, 0, 0) (approximate tilt by 
0.017° = 1') we introduce an asymmetry between the 
left and right diffraction half-plane. Thus the 'left-right' 
difference will not be a purely magnetic signal anymore 
and therefore evaluated mi,/ms ratio will contain a sys- 
tematic error. This error is surprisingly large, it even 
exceeds the asymmetry in 2BC, see Fig. |7] (middle row). 
One of its consequences is that inside the Thales circle 
the apparent mi^jms ratio in the upper half-plane is dif- 
ferent than in the lower half-plane. 

Evaluation of the 'up-down' difference leads to a result 
comparable with 2BC: there are areas of the diffraction 
plane, where the m^/ms ratio is accurately reproduced 
by the apparent mi^/mg ratio. Fig. [t] (top row). 

A quite natural idea is to perform a double difference 
procedure, i.e., to construct a 'left-right' difference map 
from the 'up-down' difference map. This procedure was 
for the first time applied in Ref. [T^l In the slightly tilted 
SBC the correction works very well in a large part of 
diffraction plane. Fig. [T] (bottom row). Note that the 
color variations are significantly reduced, meaning that 
the minims ratio can be extracted from larger part of 
the diffraction plane. Double difference procedure takes 
an advantage of the observation, that the change of the 
signal intensity caused by the tilt, is almost the same in 
the upper and lower half-plane. More detailed analysis 
of the error correction in double difference maps is given 
in Appendix [P] 

Overall, the SBC orientation provides clean magnetic 
signal and might become the prefered geometry for quan- 
titative EMCD measurements, when energy resolved 
maps of the dichroic signal are available. The small devi- 
ations from a perfect SBC orientation, can be corrected 
by performing a double difference map, thus utilizing in- 
formation from all four quadrants. 
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FIG. 8; Maps of the 'up-down' relative difference at Lz (top 
row) of 20nm thick sample of bcc-Fe as a function of the 
Laue circle center (indicated in the upper left corner of each 
panel). The corresponding rriL/ms ratio is in the middle row. 
A double difference map of the niL/ms ratio is shown in the 
bottom row. 



V. GENERAL SYSTEMATIC ROW 
CONDITIONS 

A natural generalization of the 2BC (/cc — G/2) 
and SBC {Icc = 0) systematic row geometries, where 
one studies 'up-down' difference maps, would be a gen- 
eral systematic row setup with in principle arbitrary 
Icc = 77G. We performed a set of calculations with r] 
ranging from —1/2 ('negative' 2BC) to 1 with step 0.1. 

Qualitative behavior of negative rj is similar to related 
positive ?/ values. The mirror symmetry between the ge- 
ometrical setups for ±77 lead to the following relations 
between resulting diffraction patterns: the non-magnetic 
part of the signal, which originates from real parts of 
MDFFs, does not change under mirror symmetry oper- 
ation, but the magnetic part of the signal changes (due 
to its pseudo- vector nature). Therefore, in theory, if we 
would have a separation of the diffraction pattern into 
its non-magnetic part (Re[MDFF]) and magnetic part 
(Im[MDFF]) at certain 77, we can obtain the diffraction 
pattern for —77 by mirroring both components, changing 
the sign of the magnetic component and adding them 
together. As a result, however, the total signals for ±77 
are not mirror images of each other and the difference 
is caused by magnetism. That actually provides another 
way, how the magnetic signal could be measured: first 
we measure a spectrum at tilt 7;G and detector posi- 
tion [kx, ky) and second spectrum at tilt — 77G and detec- 
tor position {—kx,ky). Their difference would be a pure 
EMCD signal. This scheme could be also extended to 
a double-difference-like procedure by measuring several 
spectra at {zLkx,iky) detector positions for both tilts. 

In Fig. [8] we show difference maps of r] — 
— O.S, 0, 0.S, 0.7 and 1. Also the apparent niL/ms maps 
are included. From these maps one can deduce that the 
'exact 2BC geometry' does not play any special role, con- 
trary to SBC where an additional mirror plane can be 
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FIG. 9: (color online) The top row shows maps of the L3 
relative difference signal of bcc-Fe in the zone axis orientation 
as a function of thickness. Bottom row shows corresponding 
mh/ms maps. 
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FIG. 10: (color online) The same as fig. [9j but with Icc = 
(100), i.e., the 2BC-like geometry close to the ZA orientation. 



used. Indeed, for 'up-down' difference maps the value of 
Icc is-within certain limits-not important. By tilting the 
beam, the area with low influence of asymmetry moves 
along with the tilt. However, the thickness dependence 
does not allow any general conclusions about optimal de- 
tector positions (see also Appendix [b| . 

In the bottom row of the Fig. [8] we show the m^/m^ 
ratio map evaluated using the double difference proce- 
dure. A somewhat surprising conclusion can be made: 
the double difference procedure seems to work reason- 
ably well also for considerable tilts. Even at 77 = 0.7 the 
double difference procedure enlarges the area with low 
influence of asymmetry. Note also the similarity between 
r] = ±0.3. Maxima at 77 = 0.3 coincide with minima at 
T] = —0.3 and vice-versa. The reason is that the asym- 
metric component of the signal changes sign, as can be 
expected from the symmetry considerations about the 
non-magnetic part discussed above. If the double differ- 
ence procedure would be applied on a sum of diffraction 
patterns at 77 = —0.3 and 77 — 0.3, the asymmetric com- 
ponent would cancel out. 



VI. ZONE AXIS CONDITIONS 

Up to now, all proposed experimental setups were de- 
vised in a way to reduce the dynamical effects as much 
as possible-in order to make the effect physically more 
transparent. Intuitively, one might expect that a large 
number of contributing MDFFs leads to concurrent phase 
shifts and thus to the reduction of the net magnetic sig- 
nal. Tilting the sample reduces the number of Bragg 
spots in the diffraction plane and thus reduces the num- 
ber of reflections contributing to the total signal. A fur- 
ther reduction of complexity of dynamical effects was ex- 
pected from the 2BC.^ 

Performing a measurement with incoming beam paral- 
lel to the zone axis leads to strong dynamical effects. A 
larger number of diffracted beams appears in the diffrac- 
tion plane. The calculation becomes considerably more 
demanding, because many points of the reciprocal lattice 
are close to the Ewald sphere and often larger number of 
Bloch coefficient products need to be taken into account. 



Still, there is no conceptual difference in the treatment 
of the zone axis case or systematic row case. 

The zone axis geometry has some advantages compared 
to 2BC or 3BC geometries. The beam tilt is zero, there- 
fore the diffraction pattern may be of higher symmetry. 
Assuming a cubic structure of the sample and (001) zone 
axis, the Icc of the incoming beam can be either zero or 
lying along any of the symmetry axes ((100) or (110)) and 
there will still be a mirror axis available in the diffrac- 
tion plane allowing to separate out the imaginary part of 
MDFFs. E.g., we can arrange a 2BC-like geometry set- 
ting the Icc to G/2 [with G = (200) or G = (110) for bcc- 
Fe] . The higher symmetry of the setup allows to take ad- 
vantage of available symmetry operations to improve the 
effective count rate and thus the SNR. When Icc = (000) 
we have four symmetry operations (90° rotations along 
the transmitted beam) connecting possible detector po- 
sitions to the symmetrically equivalent positions in the 
diffraction plane. This allows to symmetrize the signal 
and double the SNR. For the (111) zone axis the diffrac- 
tion pattern would acquire 3-fold rotational symmetry. 

We performed two sets of zone axis simulations - an ex- 
act zone axis orientation (ZA), Icc = (000) (see Fig. |9|, 
and a 2BC-like orientation (ZA-2BC) with Icc = (100) 
(see Fig. [To]), respectively. Out of these, the ZA orienta- 
tion is of higher symmetry - there are four mirror axes, on 
which the dichroic signal vanishes. All four mirror axes 
allow clean extraction of the magnetic part of the sig- 
nal (no asymmetries) , therefore it does not matter which 
one is chosen. Because the Thales circle positions for 
both G = (200) and (110) are lying on the mirror axes, 
obviously these are not suitable for detector placement .ESl 
In the ZA-2BC geometry there is only one mirror plane 
left and the Thales circle positions may be, in principle, 
used. Though they are again not optimal, as can be seen 
from Fig. [l^ 

In both cases the evaluated mL/ms ratio is smooth 
and constant throughout the diffraction plane, if one 
avoids areas with zero or low magnetic signal, where 
the error bar of 771^/7715 diverges. The signal intensity 
(and its magnetic content) naturally decreases with dis- 
tance from the transmitted beam. Due to pronounced 
dynamical effects, the magnetic signal forms a complex 
pattern with many local minima and maxima, changing 
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FIG. 11: (color online) Slightly misoriented zone axis setup 
with Icc = (0.07,0.04,0). Top row: diffraction patterns at 
Z/3 edge (arbitrary units). Note the slight differences in in- 
tensities of {110} Bragg spots, best visible at 30nm. Next 
rows show the apparent rriL/ms ratio constructed using: hor- 
izontal mirror axis (2°"* row), double difference (S''* row) and 
triple difference, i.e., using all available mirror axes (4"^ row). 
Note, how the range of the color bar can be reduced, as we 
apply more symmetrization to the data. 



sign throughout the diffraction plane. This is a potential 
disadvantage of the ZA geometry, because the magnetic 
signal is low also in the surrounding of contour lines with 
zero magnetic signal. (This also limits the size of de- 
tector aperture.) By measuring the reciprocal maps it 
is possible to identify places with strong magnetic signal 
and use these areas for m^/ms evaluation. 

In ZA-2BC geometry the magnetic signal has less sign 
changes in the diffraction plane, therefore it might seem 
to be more suitable for experiments. Though, it would 
be desirable to study the effect of small misorientation, 
particularly the influence of nonzero component of Icc 
perpendicular to G. In ZA-2BC there is no additional 
mirror axis available, which would help to correct the 
map of magnetic signal by the double difference proce- 
dure. 

In ZA orientation there are four mirror axes, which 
might be used to amplify the magnetic signal in difference 
maps in a way similar to the double difference procedure 
for SBC. We tested it on a slightly tilted ZA orientation. 
The Icc was set to (0.07, 0.04, 0) - position, which is out- 
side of any mirror axes. This introduces a slight asym- 
metry into the diffraction pattern, see Fig. |11[ Asymme- 
try seems to be rather weak, it is best seen at thickness 
30nm as an intensity variation among individual {110} 
Bragg spots. Despite the fact that it is barely visible in 
the diffraction pattern, its influence on apparent m^/m^ 
is very strong. See the map of apparent m^/ms ratio 
constructed from 'up-down' difference maps in Fig. |11[ 
It is obvious that this map would not allow to obtain 
reasonable quantitative magnetic information. A similar 



situation is obtained also for 'left-right' based maps (not 
shown) . The double difference map using both horizontal 
and vertical mirror axes considerably improves the map. 
Further improvement is achieved by exploiting also the 
diagonal mirror axes - a 'triple difference'. Nevertheless, 
there are still quite strong deviations from the correct 
mi^/nis value. As a whole, these simulations indicate 
that zone axis measurements of the rriL /ms ratio are very 
sensitive to accurate orientation of the sample. Though, 
when the orientation can be set with sub-0.1 mrad ac- 
curacy, ZA measurements will provide highly symmetric 
maps and allow an accurate 771^/7715 extraction. 



VII. CONCLUSIONS 

We explored the role of symmetry in EMCD experi- 
mental geometries, where the sample itself serves as a 
beamsplitter and the dichroic spectra are obtained as a 
difference from two detector positions. Various sources 
of asymmetries and their influences on the extracted 
nT'L/ms ratio are analysed. We have provided a survey 
of EMCD experimental geometries, stating their advan- 
tages and disadvantages. We have demonstrated, that 
the Thales circle positions are, in general, not optimal 
for detector placement. Optimal detector positions are 
sample-, thickness- and orientation-dependent. 

Our analysis has shown, that the two-beam case is, 
strictly speaking, not suitable for the sum-rules applica- 
tion due to its inherent asymmetry. Despite this asym- 
metry, the dichroic signal can be rather strong and an 
accurate m^/ms evaluation can be performed within a 
reasonable range of detector positions. This geometry 
seems to be advantageous, when it is not possible to ac- 
quire energy-resolved images of the diffraction pattern 
and only spectra at particular detector positions are mea- 
sured. We have shown that exact setting of the two-beam 
case is not crucial and deviations from two-beam case to- 
wards three-beam case do not qualitatively influence the 
tol/tos evaluation. 

We have analysed the three-beam case, which enables 
two symmetric detector positions, provided there is a 
mirror symmetry plane perpendicular to the systematic 
row index. Possible deviations from perfect three-beam- 
case geometry can be corrected by constructing a double- 
difference map. This geometry might become the method 
of choice, when energy-resolved diffraction patterns are 
available. Inf ormat ion from the whole diffraction plane 
can be utilizecP^'^l], 

Finally, as it turns out, measurement in the zone 
axis orientation is also possible. It brings an advan- 
tage of higher symmetry of the diffraction pattern, which 
may be used to improve the signal to noise ratio. On 
the other hand, pronounced dynamical diffraction effects 
cause a high sensitivity to even subtle deviations (few 
tenths of miliradians) from the exact zone axis orienta- 
tion (Icc = 0). These errors can be corrected to some 
extent by double or triple difference procedures, never- 
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theless in zone axis orientation this approach is not as 
efficient as in the three-beam case. 
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APPENDIX A: SUMMATION OVER BLOCH 
COEFFICIENTS 



Evaluation of the double differential scattering cross- 
section involves calculation of the following sum [see, e.g., 
Ref. H Eq. (24)] 



dndE 



1_ ^^u(cL_qV^ i(q_q').u 



ghg'h' 

jij'i' 



q2qi2 



where 



^ghg'h' ~ ^0 '-'g ^0 



(Al) 



(A2) 



Here C^^^ and D^'' are the Bloch coefficients for the in- 
coming and outgoing Bloch fields, g,h are Bragg reflec- 
tions and j, I are indices of the Bloch waves. The quantity 
is the MDFF divided by squares of momen- 
tum transfer vectors q, q' (Coulomb potential factors). 
For the rest of the notation we refer the reader to our 
above-mentioned paper. 

The main difficulty originates in the eight-fold summa- 
tion over Bragg reflections and Bloch-wave indices (to be 
performed Nt x Ne times, where Nt is number of sam- 
ple thicknesses and Ne is number of energy-loss steps 
considered). Since the rriL/ms ratio depends sensitively 
on accuracy of simulations, a large number of terms is 
needed to achieve good convergence of results. In our 
previous simulations a subset of Bragg reflections were 
selected by means of Wg.max and using these we selcected 
a subset of Bloch waves with certain minimum norm over 
these reflections. That was performed for both incom- 
ing and outgoig Bloch fields and all cross-terms between 
these indices were taken into the summation. This way a 
vast number of negligible terms were summed consuming 
valuable computing time. 



Here we try to pre-identify the important terms and 
evaluate only those, which are bigger than some mini- 
mum value cut-off, Pniim in this paper set to 10^^ (it 
may be increased or reduced, depending on accuracy de- 
mands). Since the thickness-dependent function Tjijui^t) 
has for given thickness t moreless constant magnitude 
and the same is true for the phase factor e^^'^~'^ their 
precise values need not to be considered in estimation of 
term magnitudes. For the MDFF, for the sake of relative 
magnitude estimation, we assume a dipole type asymp- 
totics, i.e., 

Siq,q',E) _ Aq ■ q' + B ■ (q X q') ^ qq^ 
S{qo,cio,E) ^ Aqo-q'o ^ 9o 

where A and B are system dependent model parame- 
ters and we estimate them to be of similar magnitude. 
For non-magnetic centrosymmetric samples B = 0. Here 
qo = kout - kin and qgh = qo + h - g. 

First we solve large eigenvalue problems for incoming 
and outgoing Bloch fields with Bragg-reflection cut-off 
condition Wg < 10^ That is fulfilled by 600-700 Bragg 
reflections (in bcc Fe) and the same number of Bloch 
waves. From the resulting eigenvectors (Bloch waves) 

we construct a list of Cq^*C'^'^ and Dq'^* products 
bigger than Pmin/Q, where Q is another convergence pa- 
rameter explained below. Number of elements in these 
lists varies between few thousands for systematic row ge- 
ometries up to 2-3 tens of thousands elements in zone axis 
simulations. These lists are sorted by magnitude using 
the Quicksort algorithm. 

The next step is generation of a list of quadruple prod- 
ucts of Bloch coefficients, d'^^ d':^^* / q^^,, big- 
ger than Pinin- This list is also sorted by QuickSort 
algorithm, but now by q-vectors as the primary criterion 
(i.e., grouping the terms with the same q-vectors next 
to each other) and magnitude as the second. Number 
of elements of this list is typically of the same order of 
magnitude as any of the lists of Bloch coefficient double- 
products constructed in previous step. Here it also be- 
comes clear, what is the purpose of the parameter Q: 
since the magnitude of the ratio qo/q^h can in principle 
exceed unit, it may happen that a quadruple product 
with this factor will exceed Pmin even if some of double 
products, of which it consists, are below this treshold. 
This may happen for larger scattering angles or higher 
energy losses and we adopt a very safe value Q = 10. 
It is not too costly, since the lists of double products 
are typically smaller than other lists constructed in next 
steps. 

Finally, a list of octuple products of Bloch coefficients 
Y-'!^ — — — is generated. Thanks to the way, how the 

ghg'h' i}gh9g'h' ° 

list of quadruple products is sorted, we can directly con- 
struct and output a q, q'-ordered list of octuple products 
^ghg''h" along with all jyO) (see Ref. 4), into a binary 
file for later summation. Generation of this list is linear 
in the number of generated octuple products. 
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Here we note that the factor — ^ — is the dipole-type 

asymptotic estimation of the relative magnitude of the 
MDFF with respect to the DFF, including Coulomb fac- 
tors (A more strict estimation would be possible, 
since dipole approximation underestimates the decay of 
the dipole-allowed transitions for large q vectors.) 

The list of q, q' is written to a file for an external pro- 
gram, which calculates the MDFFs. Obviously, the final 
summation is linear in the length of the list of octuple 
products, and is very fast and more accurate compared to 
a summation over all cross-terms performed before, be- 
cause with this algorithm we can afford to include much 
larger number of beams and Bloch waves. The comput- 
ing time scales moreless proportionaly to the inverse of 
^'min, which is an enormous improvement over x 
of the previous method (per thickness and energy; A^g is 
number of Bragg reflections and Nj is number of Bloch 
waves included in summation). 
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FIG. 12: Diffraction patterns (top row), relative up-down dif- 
ference maps (middle row) and apparent niL/ms ratio maps 
(bottom row) a function of convergence parameter Pmin. Cal- 
culations were performed for 20nm layer of bcc iron in two- 
beam geometry with G — (200), at SOOkeV. 



Pmin 


{N2} 


{N,} 






t 


10" 


■2 


230 


72 


330 


12 


2h 16min 


10" 


3 


1355 


750 


4710 


110 


3h 18min 


10" 


4 


5545 


5865 


56830 


705 


llh 13min 


10" 


5 


23060 


39760 


570250 


3600 


83h 35min 


10" 


6 


73591 


259739 


5398841 


16860 628h 14min 



TABLE II: Average lengths of double {N2), quadruple {A4) 
and octuple (Ns) product lists, average number of momentum 
transfer diads per energy step ( Aqq/ ) and computing times for 
maps in Fig. [12] as a function of convergence parameter Pmin. 
Times refer to a single Intel Pentium 4 Xeon processor at 
2.5GHz. 

We summarize the properties of the described algo- 
rithm in the Fig. [12] and Table [II] for five different values 
of convergence parameter Pmin- The calculation is per- 
formed on a grid of 51 x 51 outgoing beam directions 
for Ne = 10 and A^^ = 100 in the two-beam orientation. 
Note that the Nt only marginally influences the total 



computing time, since it only scales the final summation, 
which is typically the fastest part of the calculation. The 
Ne infiuences linearly the MDFF calculation and final 
summation time. Generation of octuple products and 
q, q' lists is Nt independent and only marginally Ne de- 
pendent (it only infiuences the output of q, q' list). 

Note that the diffaction pattern appears to be rea- 
sonably converged already for Pmin — 10^^. An atten- 
tive reader might spot slightly fuzziness, but neverthe- 
less, the difference between all five diffraction patterns 
are visually negligible. The relative difference map re- 
quires better convergence, at least Pmin = 10"'^, or better 
Pmin = 10~^. Note the numerical noise at Pmin = 10""^, 
particularly along the vertical line going through the 
transmitted beam. Although a slight fuzziness remains 
at Pmin = 10"'^, but the accuracy is already satisfac- 
tory. The most sensitive quantity presented here is the 
apparent m^^/mg ratio. Since the orbital momentum 
in iron is small, the nominator in the sum rule expres- 
sion is a difference of two (typically small) differences of 
spectra, integrated over L2 or L3 edge, respectivelly. It 
requires high accuracy to obtain well converged maps. 
Calculation with Pmin = lO^"* seems to produce reason- 
ably converged results. The only visible difference at 
Pmin = 10~^, when compared to Pmin = 10"'', appears 
around (— 1, ±l)G'(2oo) positions. The results of calcu- 
lation with Pmin = 10~^ are visually indistinguishable 
from Pmin — 10^^. So we can conclude that for purposes 
of this paper the convergence criterion Pmin = 10^'^ pro- 
vides highly converged data - particularly when we con- 
sider a region of sizable difference signal, where an appre- 
ciable the magnetic signal to noise ratio can be expected. 

The Table lU summarizes somewhat more technical in- 
formation about the computational complexity. We show 
the average lengths of lists of Bloch-coefhcient products, 
which scale moreless inverse-proportionally to the conver- 
gence parameter. The computing time apparently grows 
slowlier, but one should take into account a constant ev- 
erhead in these calculations connected with two diago- 
nalizations of complex matrices of dimension ca. 630 and 
the start-up and initialization costs of MDFF calculation, 
which are both repeated 51 x 51 times. These overheads 
are estimated to be roughly around 2 hours. After sub- 
traction of that a moreless linear scaling of the compu- 
tation time with inverse Pmin is observed. We note that 
the computing time is typically by factor 10-15 longer for 
zone axis calculations. 

The whole algorithm is trivial to parallelize, since each 
pixel in the diffraction plane is calculated independently. 
Therefore practically a linear scaling with number of pro- 
cessors is achieved. For applications, where only the 
diffraction pattern is needed, the method delivers accu- 
rate results at affordable computing time. 
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APPENDIX B: A NOTE ON OPTIMAL 
DETECTOR POSITIONS 

An important corollary of the sum-rules is, that the 
imaginary part of energy integrated MDFF does contain 
the magnetic information regardless of the relative ori- 
entation of q and q' momentum transfer vectors. The 
only exceptions are parallel or antiparallel orientations, 
respectively, when the imaginary part is zero (within 
dipole approximation). This implies, that the detector 
placement is by no means limited to the Thales circle 
positions (Fig. [s]). In other words, the Thales circle po- 
sitions have no special importances^. We have presented 
a set of simulations varying the detector position to scan 
a part of the diffraction plane, which lead to difference 
maps related to magnetic signal. These maps provided 
several findings: i) the maximum of the dichroic signal 
is, as anticipated, at places other than at the Thales cir- 
cle positions; ii) there can be several local maxima and 
minima in the diffraction plane; iii) the pattern of the 
dichroic signal sensitively depends on the sample thick- 
ness and orientation. 

This means, that the question of the optimal detector 
positions actually does not have a unique general answer. 
There are orientations and thicknesses, when the optimal 
detector positions are quite different from Thales circle 
positions. All this is also sample dependent. Therefore 
in actual experiments it would be advantageous to have 
access to such maps of the dichroic sig nal. Su ch experi- 
ments actually have been already done .1^21^^122] 

Another application of the above-mentioned corollary 
of sum-rules is the possibility to integrate the dichroic 
signal over some part of the diffraction plane. Because 
the prefactors in spin and orbital sum-rules are the same, 
the signal from a range of detector positions can be col- 
lected. This can be done either by numerically integrat- 
ing within the map of dichroic signaP^ or by having a 
suitable collection aperture. This observation is related 
to the question of optimizing the SNR - selection of the 
optimal size, shape and placement of the spectrometer 
entrance aperture.^i^ [Alternative approach of improving 
SNR of a measured physical quantity - mi^jras ratio - 
was developed in Ref. [121 where a statistical analysis was 
performed within mL/ms maps.] 
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FIG. 13: (color online) Mean value of the spin and orbital 
operators for unoccupied states as a function of the upper 
integration limit E. The energy range inside the blue dashed 
rectangle was considered in calculations. Ground state values 
are marked by horizontal lines. 



equally spaced energy values. This interval was chosen to 
encompass the main peak in the unbroadened spectrum 
and applies for both L2 and L3 edges. Consequences of 
this approximation are briefly discussed in this appendix. 

In Fig |13| we present the expectation values of the spin 
and orbital momentum operators evaluated as a function 
of the energy integration range, {Ep^Ep + E), within 
the d-states subspace. For infinite E these values would 
converge to ground state spin and orbital moments of 
d-states, which are represented by horizontal lines. The 
evolution of both spin and orbital momentum as a func- 
tion of energy is different, therefore also the 'apparent' 
TfiL/ms ratio will depend on the energy range. Particu- 
larly, for the range {Ep + 1.5eY , Ep + 2.0eV) the mL/ms 
ratio is 0.045 - that is roughly twice the ground state 
value. 

It is important to realize that this fact does not influ- 
ence any of the previous qualitative conclusions, which 
are system-independent. Indeed, the maps of dichroic 
signal as a function of energy show identical trends with 
varying amplitudes. 



APPENDIX C: ENERGY INTEGRATION RANGE 



APPENDIX D: DOUBLE DIFFERENCE 
CORRECTION 



Extensive calculations presented in this work would 
not be feasible without approximations, which were pre- 
sented in the description of the computational methods. 
These approximations were carefully chosen to provide 
correct unbiased results. There is one exception though, 
which is the energy integration range in m^/ms evalua- 
tion. This should be in principle infinite, {Ep, 00). Since 
the MDFF evaluation is the most time-consuming step 
of calculations, the energy integration range was reduced 
to a small interval {Ep + l.be\,Ep + 2.0eV) with 11 



The 'left-right' difference signal extracted from a per- 
fect SBC geometry was shown to be of purely magnetic 
origin. In practice, though, there is always some devi- 
ation from perfect SBC towards G or — G. [The beam 
tilt leading to a systematic row setup is of much larger 
magnitude (^ 0.1 rad) than tilts discussed here (sub- 
miliradian). Therefore deviations along ky do not play 
role-they are simply absorbed into the beam tilt-and one 
only has to deal with deviations along (the x axis is 
along the systematic row).] It was demonstrated in Fig.m 
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FIG. 14: (color online) Map of the 2ARc/Aim for a slightly 
tilted SBC, Icc — (0.05,0,0). Top row: 'up-down' difference 
(horizontal mirror axis), middle row: 'left-right' difference 
(vertical mirror axis), bottom row: double difference (both 
mirror axes). 



that double difference correction significantly extends the 
area, where an accurate m^/ms ratio extraction is pos- 
sible. Here we will briefly discuss, why this correction 
works. 

By slightly tilting the beam towards positive k^, the 
intensity throughout the diffraction pattern increases in 
the right half-plane at all energy levels. Particularly, this 
happens for the nonmagnetic component of the inten- 
sity. We assume that the increase is similar in upper 
and lower half-planes, i.e., Al^k^, ky) « A/(fc^,— fcy), 
while AI{—kx,ky) is of different sign. Then the left- 
right difference will contain an error proportional to 
AI(kx, ky) — A/(— fca;, fcj,), which increases the magni- 
tude of magnetic signal on one edge and decreases it 
on the other edge (assuming that they are of different 
sign). This biases the observed rriL/ms value. Addition- 
ally, this error influences the upper and lower half-plane 
in opposite ways and thus the extracted mL/ms ratio 
differs in upper and lower half-planes. 

If we perform a 'left-right' difference, the nonmag- 
netic components from the exact SBC orientation can- 
cel out and only the "error" superimposed on magnetic 
signal remains. Subsequent 'up-down' difference en- 
hances the magnetic signal and, simultaneously, reduces 
the error thanks to its approximate up/down symmetry 
mentioned above: AI{kx,ky) ~ AI{kx,—ky). Remain- 
ing error, [AI{kx,ky) - AI{kx,~ky)] - [AI{-kx,ky) - 
AI(—kx,—ky)], is considerably smaller than asymmetry 



in pure 'up-down' or 'left-right' differences, respectively. 

To verify these assumptions quantitatively we calcu- 
lated separately the magnetic and nonmagnetic contribu- 
tion to the scattering cross-section for the slightly misori- 
ented SBC with Icc = (0.05,0,0). This amounts to two 
separate calculations considering either only real parts 
(non-magnetic component) or imaginary parts (magnetic 
contribution) of MDFFs, respectively. 

The difference maps of the non-magnetic part of the 
signal (Aro; not shown) are directly related to the error 
introduced by asymmetries. Very instructive is to show a 
ratio Arc/Ai,„, i.e., the difference of non-magnetic signal 
divided by the corresponding difference of the magnetic 
signal. Such maps describe the relative enhancement of 
the apparent magnetic signal due to non-magnetic contri- 
bution. Maps of ARo/Aiin evaluated at the L3 edge are 
shown in Fig.[T4]for 'left-right', 'up-down' and double dif- 
ference, respectively. Maps at L2 edge (not shown) have 
the same trends, but with opposite sign and reduced mag- 
nitude. Note that values in the double difference map are 
much lower than in the other two maps, demonstrating 
the effectiveness of the the double difference procedure. 

Note that these maps are strikingly similar to Fig.[7j It 
should not be surprising though, since they quantify the 
error caused by asymmetry. The error introduced into 
ttT'l/ttT'S evaluation due to asymmetry can be expressed 
more precisely in terms of the following equations. Let's 
assume that the errors A^e at L3 are twice as large com- 
pared to errors at L2 edges. This is approximately the 
case for late Sd metals, where the branching ratio is close 
to 2. Then the apparent m^/ms ratio is 

2 (Al3 + 2ARe) + (Al, + Aro) 



S (Al3 + 2ARe) 

2 Al,+Al 

S Ar, - 



+ SAro 



Arc) 



ms 



-2Ai, 

2ARe 



2A, 



ms 



2A 



Ro 



3Ai„ 



stands for Aim at particular edge. In the 



last step we assumed that A^^ 



This 



equation expresses the quantitative influence of the error 
caused by asymmetry in terms of the ratio ARe/Ajm, and 
explains the similarity of Fig. 14 and Fig.jTj We note that 
the error is additive. 
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Since crystal is tilted the outgoing beam directions at de- 
tector positions (A) and (B) have different components 
perpendicular to the sample surface and thus the elonga- 
tions 7'-'-' of Bloch-wave vectors k'-*' will differ slightly, see 
Eq. (7) in Ref. |^ This causes a weak asymmetry in scat- 
tering cross-sections even within systematic row approxi- 
mation. Its influence depends on sample thickness and can 
be seen in 2BC maps of the rriL/ms ratio, Fig. [l] 
Within the 2BC it was tried to work with a simple ana- 
lytic model to describe the inelastic scattering using the 
transmitted beam and one Bragg scattered beam only. As 
it turned out, even in 2BC setup it is important to include 
other beams on the systematic row for quantitative EMCD 
description. And for higher thicknesses we need to abandon 
the systematic row approximation completely. 
This information may also be important for experiments 
with small tilt angles, because then the experimental setup 
is rather close to ZA orientation and thus only a weak 
magnetic signal is expected. Moreover it would be heavily 
influenced by dynamical effects. 

This relates to experiments described here, where the crys- 
tal itself is used as a beamsplitter. 



